Scaling properties of the projected SO {5) model in three dimensions 
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We study the scaling properties of the quantum "projected" 50(5) model in three dimensions by 
means of a highly accurate Quantum-Monte-Carlo analysis. Within the parameter regime studied 
(temperature and system size), we show that the scaling behavior is consistent with a 50(5)- 
symmetric critical behavior in a relative extendent transient regime. This holds both when the 
symmetry breaking is caused by quantum fluctuations only as well as when also the static (mean- 
field) symmetry is moderately broken. We argue that possible departure away from the 50(5) 
-symmetric scaling occurs only in an extremely narrow parameter regime, which is inaccessible both 
experimentally and numerically. 

PACS numbers: 74.20.-z, 74.25.Dw, 11.30.Ly, 02.70.Uu 



INTRODUCTION 

A common feature of the phase diagram of most higli- 
Tc superconductors (HTSC) is the close proximity of the 
superconducting (SC) and the antiferromagnetic (AF) 
phases. The 5*0(5) theory of high-Tc SC describes the 
transition between these two phases by an effective quan- 
tum non-linear a model with approximate 5*0(5) symme- 
try, which unifies the two order parameters Several 
microscopic 50(5) -symmetric models have been pro- 
posed which succesfullv describe many features of the 
cuprate physics [2, y, Ij; 1^ IS| ■ A breakthrough has been 
recently achieved by a model that reconciles the 50(5) 
symmetry with the physics of the Hubbard gap 0, Mi- 
In the so-called "projected" 5*0(5) (p5*0(5) ) model, the 
Gutzwiller projection is implemented exactly. Although 
the exact symmetry is destroyed at the microscopic level 
due to quantum fluctuations originating from the pro- 
jection, it would be interesting to investigate whether 
an effective 5*0(5) symmetry still controls the asymp- 
totic behavior of the system at long- wavelength, i. e., 
50(5) symmetry is asymptotically restored. This has 
been shown to be the case, for example, for a two-leg lad- 
der system. This behavior is physically similar to the 
case of interacting Fermi liquids in which the character- 
istic symmetries of the non-interacting system - namely 
conservation of the particle number for each momentum 
k of the Fermi surface - are destroyed by the interac- 
tion, but are asymptotically restored for low-energy ex- 
citations close to the Fermi surface. In contrast to the 
ladder, for a higher-dimensional system, the candidate 
for symmetry restoration should be sought at finite tem- 
peratures and is provided by the multicritical point where 
the AF and SC critical lines meet. 

The problem of symmetry restoration at this multi- 
critical point has been addressed by Hu in Ref. via a 
Monte Carlo (MC) calculation of a classical 5*0(5) model 
in which an additional quartic anisotropy term has been 



included. Classical MC are by orders of magnitude eas- 
ier to perform and less resource demanding than QMC 
simulations, hence very large system sizes can be simu- 
lated and highly accurate data are obtained. For this rea- 
son Hu could carry out a detailed analysis of the AF-SC 
phase diagram and of the critical behavior. In particu- 
lar, he could show that the SC helicity modulus follows a 
scaling law determined by the corresponding 50(5) crit- 
ical exponent z^5 and by the crossover exponent </>. For 
these exponents, Hu could find a good agreement with 
the results of the e-expansion. On the other hand it was 
pointed out by Aharony |^ via a rigorous argument that 
the decoupled fixed point is stable, and he further con- 
cluded that neither the biconical nor the 50(5) fixed 
points are stable. However, he also commented that the 
unstable flow is extremely slow for the 50(5) case due 
to the small crossover exponent. 

In this paper, we address this issue by a detailed nu- 
merical analysis of the phase diagram and of the crit- 
ical properties of the three-dimensional p50(5) model. 
Our strategy is, thus, to start from a model in which 
50(5) symmetry is realistically broken and to investigate 
to what extent this symmetry is restored at the multi- 
critical point. Our main results are twofold. Within the 
system sizes and temperature ranges we can achieve, the 
scaling behavior is consistent with an 5*0(5) critical be- 
havior. On the other hand, since the 50(5) fixed point 
could be ultimately unstable 0, we make an analyti- 
cal estimate - based on the e-expansion - of the size of 
the critical region in which deviations from the 5*0(5) 
behavior should be observed. Due to the small crossover 
exponent it turns out that the unstable flow only takes 
effect when the reduced temperature measured from the 
bi-critical point is very small. Therefore, the possible un- 
stable effect can neither be observed experimentally nor 
numerically. For all practical purposes, the multicritical 
point is dominated by the initial flow towards the 50(5)- 
symmetric behavior[l(J . Again, this situation is very sim- 
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ilar to the case of Fermi liquids. Due to the well-known 
Kohn-Luttinger effect j the Fermi-liquid fixed point is 
always unstable towards a SC state. However, this effect 
is experimentally irrelevant for most metals since it only 
works at exponentially low temperatures. 



THE MODEL 

In the pSO{5) model each coarse-grained lattice site 
represents a plaquette of the original lattice model, and 
the lowest energy state on the plaquette is a spin sin- 
glet state at half-filling. There arc four types of excita- 
tions, namely, three magnon modes and a hole-pair mode. 
Their dynamics are described by the following Hamilto- 
nian: 

3:,a = 2,3,4 X 

- Js ^ naix)naix')-Jc^{tl{x)thix') + h.C.), 
<xx'>^a=2,3,4 <xx'> 

Here t^^}^2 3 4 anihilates (creates) a triplet state, ani- 
hilates (creates) a hole pair state and Ua = {ta + ta)/V2 
are the three components of the Neel order parameter. 
As and Ac ~ U are the energies to create a magnon and 
a hole-pair excitation, respectively, at vanishing chemi- 
cal potential /x = 0. This model can also be effectively 
obtained by a coarse-grained reduction of more common 
models such as t — J or Hubbard ^3 ■ In order to study 
the effect of symmetry breaking we consider different 
situations associated with different sets of parameters. 
First, we consider the case where Js = J = (our zero 
of the chemical potential is such that A^ — Ac). It has 
been shown that this model has a static S'0(5) sym- 
metry at the mean-field level and that the symmetry is 
only broken by quantum fluctuations ^3 ■ Since we want 
to carry out our analysis also for a more realistic model in 
which also the static S'0(5) symmetry is broken, we also 
consider a system with a different ratio Jg/ Jc- In partic- 
ular, one would like to reproduce the order of mangitude 
of TcjTjq observerved in the cuprates, where Tc (TV) 
denominates the SC critical temperature (Neel Temper- 
ature). However, this behavior is obtained for Js/ Jc ^ 2, 
for which the numerical simulation is rather unstable, 
making it impossible to determine the critical exponents 
with sufficient accuracy. For this reason, we choose a 
value of the parameter "in between" (Jc = Jh = J), for 
which also the static 5*0(5) symmetry is broken. 

The phase diagram of this model in two dimensions 
has been analyzed in detail by a numerical Quantum- 
Monte-Carlo approach in Ref. 0|. In particular, the 
model has been shown to provide a semiquantitative de- 
scription of many properties of the HTSC in a consistent 
way. In Ref. the SC transition has been identified 



as a Kosterlitz-Thouless phase in which the SC corre- 
lations decay algebraically. Unfortunately, there is no 
such transition for the AF phase in two dimensions, as 
all AF correlations decay exponentially at finite temper- 
atures. Therefore, in order to analyze the multicritical 
point where the AF and SC critical lines meet, it is nec- 
essary to work in three dimensions, which is what we 
investigate in the present paper. 

RESULTS 

Case Js = Jc/2 

We start by presenting the phase diagram of the 3D 
pSO{5) model for the "symmetric" case Js = Jc/2. Fig- 
ure n shows an AF and a SC phase extending to fi- 
nite temperatures as expected. Furthermore, the two 
phase transition lines merge into a multicritical point 
(at Tb = 0.960 ± 0.005 and Hb = -0.098 ± 0.001). The 
line of equal correlation decay of hole-pairs and triplet 
bosons also merges into this multicritical point P. Un- 
like the corresponding phase in the classical model, the 
SC phase extends only over a finite /i range; this is due 
to the hardcore constraint of the hole-pair bosons and 
agrees with experimentally determined phase diagrams 
of the cuprates. In this sense, the quantum mechanical 
pS'0(5) model is more physical than the classical SO{5) 
model. However, in real cuprates the ratio between the 
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FIG. 1: Phase diagram r(^) of the three-dimensional pSO{5) 
model with J = Ja = Jc/2 and As = Ac = J. and Nt are, 
respectively, the hole-pair and the magnon-dominated regions 
of the disordered phase. The separation line between Nh and 
Nt is the line of equal spatial correlation decay of hole-pairs 
and bosons. The inset shows a detailed view of the region 
near the multicritical point. 

maximum SC temperature Tc and Neel temperature T/v 
is about 0.17 to 0.25, whereas in the p50(5) model we 
obtain the values TJJ = 1.465 ± 0.008 at Hopt/J ~ 1-7 
and Tn/J= 1.29 ± 0.01, hence Tc is slightly larger than 
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T/v. In order to obtain realistic values for the transition 
temperatures, it is necessary to relax the static 5*0(5) 
condition and take a smaller value for the ratio Jc/{2Js), 
which breaks SO{5) symmetry even on a mean field level. 
The phase diagram with Jc/{2Js) = 0.225 is plotted in 
Fig. 121 As one can see, this gives a more realistic ra- 
tio of T/v/Tc « 0.2. However, it should be pointed out 
that the numerical efi^ort to treat such different values of 
J is order of magnitudes larger than considering Jc and 
Js of the same order of magnitude, as we have done in 
Fig. ^ Therefore, we will also consider a system with 
Jc = Js = 1 for which also the static SO{5) symmetry is 
broken. For the same reason, we neglect here the c-axis 
anisotropy and consider an isotropic 3D model. 
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FIG. 2: Phase diagram for Jc/{2Js) = 0.225 as function of 
the hole doping S. 

We first carry out an analysis of the critical properties 
for Jc/ (2 Jg) = 1 A closer look to the phase transition line 
between the points S and P reveals (inset of Fig.^ that 
this line is not vertical as in the classical SO{5) model but 
slightly inclined. This indicates that a finite latent heat 
is connected with the AF-SC phase transition. Moreover, 
this means that in contrast to the classical model, ^ is 
not a scaling variable for the bicritical point P. 

Scaling analysis 

We now perform a scaling analysis similar to the one 
performed by Hu in a classical 5*0(5) system. The 
most important outcome of this analysis will be the 
strong numerical indication that in a large region around 
the multicritical point the full 50(5) symmetry is ap- 
proximately restored. This is non trivial for a system 
whose 50(5) -symmetry has manifestly been broken by 
projecting out all doubly-occupied states. First we want 
to determine the form of the T/v(/x) and Tc{n) curves in 
the vicinity of the bicritical point. For crossover behavior 
with an exponent > 1 one would generally expect that 



the two curves merge tangentially into the first-order line. 
However, this holds for the scaling variables, therefore, 
one should first perform a transformation from the old /U 
axis to a new fj,' axis defined by fJ-'{T) = fi~ (T ^Tf,)/m , 
where m « 0.11 is the slope of the first order line below 
Tb. 

After this transformation, the transition curves Tj^i^fjl) 
and Tc{n') are quite well described by the crossover be- 
havior (we now drop the prime for convenience) 

^ - 1 = 52 • (A* - A^.)^/^ 
and ZM^_i = i33.(^,_^)i/0 (2) 

The fit to this behavior is shown in more detail in Fig. El 
However, the value of (j) we obtain {(j) « 2.35) is consider- 
ably larger than the value expected form the e-expansion. 
It should be noted that the above determination of (p is 
not very accurate: the data points in Fig. (21 are the re- 
sult of a delicate finite-size scaling for lattice sites up to 
18'^, followed by the transformation from n to n' which 
again increases the numerical error bars. For this reason 
it cannot be excluded that the difference in the (p values 
is mainly due to statistical and finite-size scaling errors. 
In fact, a more accurate evaluation of tj) will be provided 
below. 

On the SC side, the finite-size scaling carried out in 
order to extract the order parameter and the transition 
temperature turns out to be quite reliable. On the other 
hand, on the AF side, the fluctuations in the particle 
numbers of the three triplet bosons slightly increase the 
statistical errors of the SSE results and make the finite- 
size scaling more difficult. 

The critical exponents for the onset of AF and SC order 
as a function of temperature for various chemical poten- 
tials can be extracted from Fig. [31 Far into the SC range, 
at /i=1.5, we find for the SC helicity modulus [l5| 

T o^{l -T/TcY with = 0.66 ±0.02, 

which matches very well the values obtained by the e- 
expansion and by numerical analyses of a 3D XY model. 
On the AF side, error bars are larger, as discussed above. 
We obtain for the AF order parameter 

Caf{oo) cx (1 - T/Tcf^ with /33 = 0.35 ± 0.03, 

for /i = —2.25, also in accordance with the value expected 
for a 3D classical Heisenberg model. 

In order to determine v and cj) more accurately in the 
crossover regime, we use two expressions derived from 
the scaling behavior (cf. Rcf. Q). 

T(r,, y.)/T{Tb, m") = ((m - M6)/(m" - t^b)y'"^ . (3) 
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FIG. 3: Plot of the AF (left) and SC (right) critical lines in the vicinity of the multicritical point. 
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where /ii, /i'j^, /i2, and /ij are related by {^i — Hhj/il^'i — 

= (^^2 - A*b)/ (M2 - ^^h) > 0. 
The result is shown in Fig. ^ we obtain the ratio 

1/5/'/' ==0.52 ±0.01, 

which is in excellent accordance with the results of the 
e-expansion and other numerical analyses 0. </> is then 
obtained by using |0J . We have applied Q onto 9 
different combinations of (jii, = /i2,Ai2) values with 
/ii/^'i =M2//i2 = 0.5. The result is 

(j) = 1.43 ±0.05, 

which is again in good agreement with the e-expansion for 
a 5*0(5) bicritical point and with the results of Ref. Q. 
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FIG. 4: Helicity T as a function of the chemical potential 
fi at T — Tt. From this function, the value of i^s/i^ can be 
extracted via equation 



as a surprise, since the 5*0(5) symmetry is only broken by 
quantum fluctuations for the parameter we have taken. 
The question we want to adress now is wether 5*0(5) sym- 
metry is also asymptotically restored for a more realistic 
set of parameters for which the static 50(5) symmetry 
is broken as well. As already mentioned above, the case, 
where the phase diagram of the cuprates is qualitatively 
well reproduced (Jc/(2Js) — 0.225, see Fig. |5J|, is too 
difhcult to address numerically, so that the critical ex- 
ponents cannot be determined with sufficient precision 
in this case. Therefore, we repeat our analysis for the 
model in an intermediate regime (Jc = J/i), which is not 
so realistic but for which the static 50(5) symmetry is 
broken as well. One could hope that if 50(5) symmetry 
is restored for here, then it might be also restored for 
the case Jc/ (2 J^) = 0.225, although one may expect that 
the asymptotic region in which this occurs will be less 
extended. We stress again the fact that eventually one 
should expect the system to flow away from the 50(5) 
fixed point, although in a very small critical region 0. 

The phase diagram for Jc ~ Jh is presented in Fig. El 
and a detailed view of the region close to the bicritical 
point is plotted in Fig. |H1 Here, the points in the plots 
were obtained by a finite-size scaling with lattices up 
to 5032 (18^) sites. In some cases, we were able to 
simulate lattices up to 10648 (22'^) sites. An example of 
the finite-size scaling is shown in Fig. [7| Our analysis 
yields = 0.682 ± 0.005 and /i;, = 0.548 ± 0.0005. Here 
the line of equal correlation decay is vertical within the 
error bars, so the transformation from ^ to /.t' is not 
necessary and the error bars are not increased by the 
transformation. This allows to determine the critical ex- 
ponents by fitting the data points visible in Fig. Elagainst 

r(^i) = Tfa * (l + {B2 + * Sign[^ib -n])*\x-iib\ 
We obtain: 



Case Js — Jc 

This agreement between the critical exponents ob- 
tained in the previsous section may not come completely 



B2 

Tb 



0.47 ±0.07, 
0.11 ±0.04,, 
1.49 ±0.18, 
0.683 ±0.004, 



(5) 
(6) 
(7) 
(8) 
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= 1.67, ±0.36 

B3 



(9) 



Since points further away from the bicritical point are 
expected to show a larger deviation from the bicritical 
behavior, we also performed a weighted fit, which takes 
this fact into account. Here, data points closer to the 
bicritical point are weighted more than the ones further 
away. Specifically, in both the SC and the AF phase, the 
point closest to the bicritical point is weighted six times 
the one with the largest distance to the bicritical point. 
The second closest is weighted 5 times and so on. The 
results are, within the error bars, quite similiar to the 
ones obtained without this different weighting procedure: 



B2 
B3 

Tb 
B2 
B3 



0.46 ±0.05, 
0.11 ±0.03, 
1.53 ±0.12 
0.682 ±0.003 

1.61 ±0.23 



(10) 

(11) 
(12) 

(13) 
(14) 



The agreement between Eas. l5l9l and Eas. lllini suggcsts 

that the data we have considered are still controlled by 
the bicritical point. 

In order to test whether alternativly proposed fixed 
points may be excluded, we carried out a least-square fit 
of our data to the decoupled fixpoint case {(/> = 1,-82, B3 
and Tf, arbitrary). This is shown in Fig. |H1 (dashed- 
dotted line). As one can see form the curve, our data 
doe not support the assumption that the pS'0(5) model 
has a decoupled fixpoint. Since wc know that ultimately 
the decoupled fixed point should be the stable one [9|, 
this means that our systems is still quite far away from 
the ultimate asymptotic region. 
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FIG. 5: Phase diagram as a function of the chemical potential 
for Jc = Jj, = 1, the lines are guides to the eyes. 
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FIG. 6: Detailed view of the phase diagram as a function 
fi {Jc = Jh ~ 1), the two lines have been obtained by fits 

to T(^) =Tb*(l + {B2 + B3* Sign[^ib - H)* I M - Mt]^) • 

The continuous (dashed) line is the 'normal' ('weighted') fit. 
The decoupled fixpoint case is plotted as a dashed-dotted line. 
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FIG. 7: Finite size scaling of the antiferromagnetic order pa- 
rameter for /i = 0.5, the temperatures cover 0.72J (lozenge), 
0.73J (star), 0.735 (square), 0.74J (triangle) and 0.75 (cirle). 
The lattice size was varied from 216 (8^) upto 10648 (22^) 
sites, scanning all cubes with even edge length. 



DISCUSSION AND CONCLUSIONS 

Altogether, the scaling analysis of the 3D pS'0(5) 
model has produced a crossover exponent which matches 
quite well with the corresponding value obtained from 
a classical SO{5) model and from the e-expansion. This 
gives convincing evidence that the static correlation func- 
tions at the pS'0(5) multicritical point is controlled by 
a fully SO{5) symmetric point, at least in a large tran- 
sient region. However, one should point out that within 
the statistical and finite-size error, as well as within the 
error due to the extrapolation of the e-expansion value 
to e = 1 one cannot exclude that the actual fixed point 
one approaches is the biconical one, which has very sim- 
ilar exponents to the isotropic 5*0(5) one. On the other 
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hand, the biconical fixed point should be accompanied 
by a AF+SC coexistence region (as a function of chem- 
ical potential), which we do not observe. As discussed 
above we can certainly exclude in this transient region 
the decoupled fixed point for which 0=10. Of course, 
our limited system sizes cannot tell which fixed point 
would be ultimately stable in the deep asymptotic region. 
Here, Aharony's exact statement shows that the decou- 
pled fixed point should be ultimately the stable one in 
the deep asymptotic region j^. 

We argue that the resolution between this exact re- 
sult and the numerically observed SO{5) critical behav- 
ior lies in the size of the critical region 0. We now 
give an estimate, based on e expansion, for the scale at 
which the instability of the SO{5) fixed point could be 
detectable. This estimate holds for the case in which one 
has a "static" SO {5) symmetry at the mean-field level. 
The symmetry-breaking effects due to quantum fluctua- 
tions have been estimated in Ref. and are given by 
Eq. (36) there. By replacing the initial conditions for the 
bare couplings in terms of the microscopic parameters of 
the Hamiltonian (cf. Eq. 26 of Ref. [iSj), and project- 
ing along the different scaling variables around the 5*0(5) 
fixed point, one obtains a quite small projection along the 
variable that scales away from the fixed point. Combined 
with the fact that the exponent for this scaling variables 
is quite small (A = 1/13 at the lowest-order in the e ex- 
pansion, although more accurate estimates 0, 0, 0| 
give a somewhat larger value of A « 0.3), we obtain an 
estimate for the scaling region in which the 50(5) fixed 
point is replaced by another - e.g. the biconical or the 
decoupled - fixed point at t={Tb- T) /T^ - IQ-i" if one 
takes the 0(e) result for the exponent. Notice that tak- 
ing the result of Ref. 0| for the exponent, one obtains a 
quite larger value t ~ 2.10"'^. However, since the multi- 
critical temperatures of relevant materials (organic con- 
ductors, and, more recently. 1^-6020^306.35) are around 
10 if, the critical region is still basically unaccessible ex- 
perimentally as well as with our quantum simulation. On 
the other hand, the other scaling variables, although be- 
ing initially of the order of 1, rapidly scale to zero due 
to the large, negative, exponents. Therefore, the 50(5) 
regime starts to become important as soon as the AF 
and SC correlation lengths become large and continues 
to affect the scaling behavior of the system basically in 
the whole accessible region. 

In conclusion, we have shown numerically that the 
pS'0(5) model, which combines the idea of 50(5) sym- 
metry with a realistic treatment of the Hubbard gap, re- 
produces the sahent features of the cuprate's phase dia- 
gram. Furthermore, this model is controlled by a 50(5) 
symmetric bicritical point, at least within a large tran- 
sient region. We have shown that this also holds for a 
case in which the static 50(5) symmetry is broken. Pos- 
sible flow away from the symmetric fix point occurs only 
within an extremely narrow region in reduced tempera- 



ture, making it impossible to observe both experimen- 
tally and numerically. We would like to point out that 
this situation is very similar to many other examples in 
condensed-matter physics. The ubiquitous Fcrmi-liquid 
fix point is strictly speaking always unstable because of 
the Kohn-Luttinger effect |ll|. But for most metals this 
instability occurs only at extremely low temperatures, 
and is practically irrelevant. Another example is the "or- 
dinary" superconductor to normal-state transition at Tc. 
Strictly speaking, coupling to the fiuctuating electromag- 
netic field renders this fix point unstable[l9j. However, 
this effect has never been observed experimentally, since 
the associated critical region is too small. Therefore, ir- 
respective of the question of ultimate stability, we argue 
that the 50(5) fix point is a robust one in a similar sense, 
and it controls the physics near the AF and SC transi- 
tions. 
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